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ABSTRACT 

The differential equations and boundary conditions describing the 
periodic-flow heat exchanger (frequently referred to as rotary regen- 
erator) are sufficiently complicated to preclude complete analytical 
solution. Because of the interest in this type of heat exchanger con- 
siderable effort has been expended to obtain useful solutions to the 
analytical problem of predicting its performance. Several approximate 
solutions have been proposed but each of these is limited due to the 
restricted ranges of parameters in which they apply. 

A numerical finite-difference method of calculating the spatial 
temperature distribution in the periodic-flow type heat exchanger is 
presented by considering the metal "stream" in cross-flow with each of 
the fluid streams as two separate but dependent exchangers. In the 
development no assumptions are made which would restrict the range of 
parameters in which the analysis would be applicable. The exchanger 
effectiveness has been evaluated over the following ranges of dimen- 
sionless parameters: 

1.0 2 Cain/Cnax % 0.10 
1.0 < Cr/Cmin < co 
1.0 < NW < 10 
1.0 > (hA)' > 0.25 

The National Cash Register 102A general purpose digital computer 
was used to accommodate the large number of subdivisions necessary for 
accuracy and an extrapolation of data to zero element area (infinite 
subdivision) was used to obtain values of effectiveness good to four 
Significant figures. 

The author expresses his appreciation to Assistant Professor C.P. 


Howard for his direction and encouragement in this work. 
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NOMENCLATURE 


SYMBOLS 
A = heat transfer area on side designated by subscript, sq ft 


C 


1 


heat capacity rate (Wc) of fluids or rotor matrix according to 
subscript, Btu/hr, deg F 
с з specific heat of fluids (constant pressure) or rotor matr@x 
material depending on subscript, Btu/lb, deg F 
h = unit conductance for thermal convection heat transfer, 
Btu/hr, sq ft, deg F 
k = unit thermal conductivity, Btu/hr, sq ft, deg F/ft 
K = a constant, defined as used and distinguished by numerical 
subscript 
N = number of subdivisions of fluid or matrix solid flow "stream" 
according to subscript 
Q = heat transfer rate, Btu/hr 
T = temperature of fluid or matrix, inlet or outlet, depending D 
subscripts, deg F 
W = mass flow rate of fluid (lb/hr) or matrix (rev/hr)(mass) 
according to subscript 
DIMENSIONLESS PARAMETERS 
E = exchanger heat transfer effectiveness; ratio of actual to 
thermodynamically-limited maximum possible heat transfer rate 
Cmin/Cmax = capacity rate ratio of fluid flow streams 
Cr/Cmin = capacity rate ratio of rotor matrix to minimum fluid 
capacity rate 
(hA)' = (hA)n/(hA),, geometric symmetry and heat transfer balance 


of exchanger 





NTU = (hA)/C，number of transfer units on side designated by 


subscript 


— 


/ 


----- [-------------- | - — 
NIUo = ος p: f / (5А), / = ντα. | η ) 


over-all number of transfer units 
SUBSCRIPTS 


avg - average 


r = rotor matrix 
1 = inlet 
ο = outlet, or over-all in case of NTU, 


min = minimum magnitude 
max = maximum magnitude 

n = value on side of Chin 
x = value on side of Cmax 
ce = cold side 

h = hot side 


co = subscript on E to indicate value when Cr/Cnin = X 


vi 





1. Introduction 

In the past decade considerable interest has been shown in the ap- 
plication of the periodic-flow heat exchanger (rotary regenerator) to the 
Gas-turbine cycle for exhaust-gas thermal-energy regeneration. Regen- 
eration by whatever method applied to the gas-turbine cycle obtains, among 
other improvements, higher thermal efficiency, lower specific fuel con- 
sumption, improved part load performance and lower optimum pressure ratio. 
Particular advantages to be expected from the application of the periodic- 
flow heat exchanger to the gas-turbine process are discussed by Harper 
and Rohsenow (1)* and Coppage and London (2), among others. Principal 
among these is the possibility of a more compact unit than is possible 
with the conventional stationary surface types. This is due to the fact 
that the periodic-flow heat exchanger matrix can be finely divided, with- 
out structural difficulty, resulting in an appreciably reduced heat ex- 
changer volume relative to a stationary type for a given effectiveness 
and pressure drop. The major disadvantage of the rotating periodic-flow 
type is the sealing difficulties encountered and the attendant leakage 
involved when an appreciable pressure ratio is used. This has been re- 
duced to acceptable limits in at least two current development programs 
(3,4) and further refinements should follow. 

Because of the interest in this type of heat exchanger there has 
been considerable effort expended to obtain useful solutions to the 
analytical problem of predicting its performance. The differential 
equations and boundary conditions describing the system are sufficiently 


complicated to preclude complete analytical solution. An analytical 





Numbers in parentheses refer to the Bibliography at the end of the thesis. 





solution is available for the special case of infinite matrix rotative 
speed where the expression for effectiveness takes the same form as that 


of a counter-flow direct type (2) and is given by 
“NTU, СІ Cm» /Cmax ) 


ATP . 
E E | A Cnn Cn ax e” MTUo (17 (min /Cmax) [1] 


and for the further condition of Cnin/Cmax = unity, reduces to 
` 





Ε = NTL。 lta] 
= | ΕΝΤίο 


Consequently, several approximate methods of solution have been 

proposed. Perhaps the oldest is that of Hausen (6), described and elab- 
orated by Johnson (8,9), which begins with the known solution for the 
initial cooling of a matrix element, initially at uniform temperature, 
by the passage of a cold fluid. By dividing the regenerator tube into a 
large number of strips of equal length and applying this elemental solu- 
tion to each strip in turn, a family of simultaneous linear algebraic 
equations may be derived whose solution gives the desired temperature 
variation in the periodic case. 

A second method, due to Iliffe (10), uses the solution to the prob- 
lem of the initial cooling of a matrix element at arbitrary initial temp- 
erature distribution, derived by Nusselt. Using this elemental solution, 
the method of Iliffe leads to two simultaneous integral equations which 
require numerical methods for their solution. In a third method, Saun- 
ders and Smoleniec (7) replace the governing partial differential equations 
by their equivalent finite difference equations and effect an approximate 
solution by relaxational techniques. In a fourth approach, Tipler (16) 
replaces the thermal problem by its electrical analogue and observes the 


solution on an oscilloscope. In yet a fifth approach, the perturbation 





method used by Schultz (8) has been extended by Jones and Fax (5) lead- 
ing to algebraic statements for the departure of the regenerator effec- 
tiveness from its asymptotic value, given by equation D] , for several 
types of asymmetry. Each of these solutions is limited in one respect 
or another, primarily due to the restricted ranges of parameters in which 
they apply. 

It is the purpose of this thesis to present the results of a finite- 
difference numerical analysis which was carried out with the aid of a 
digitial computer such that there were no limitations imposed on the four 


independent dimensionless parameters selected. 





2. Method 


The conventional idealizations and boundary conditions assumed in 


the usual derivation of the governing differential equations (2) are as 


follows: 


a. 


The thermal conductivity of the matrix is zero in the direction 
of fluid flow and matrix metal "flow". It is infinite in the 
other direction normal to the fluid flow. 

The specific heats of the two fluids and matrix material are 
constant with temperature. 

No leakage of the fluids occurs either due to direct leakage 
or carry over, and each fluid flow is unmixed. 

The convective conductances between the fluids and the matrix 
are constant with flow length. 

The fluids pass in counterflow directions. 

Entering-fluid temperatures are uniform over the flow inlet 
cross section and constant with time. 

Regular periodic conditions are established for all matrix 


elements, i.e., steady state condition. 


With these idealizations a typical element, Fig. 1, from each fluid 


side of the exchanger can be represented schematically as shown in Fig. 


2. Although Fig. 1 is a disk (axial-flow) type exchanger, typical el- 


ements of a drum (radial-flow) type can be selected, with the same ass- 


umptions, which can be represented schematically in the same manner as 


Fig. 2. 


It is apparent from Fig. 2 that each of these elements can be regard- 


ed as a crossflow exchanger with a gas stream and a metal "stream". This 


extension of the method of calculating crossflow exchangers to the rotary 
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regenerator was suggested бу С. Μ. Dusinberre in his discussion of ref- 
erence (2). It should be noted that the Tp; and Ty, on the two elements 
are not the same, but to avoid further subscript notation the side in 
question will be specified when it is not obvious. However, the Ty. of 
a particular element is the Ty; of the next adjacent element in the 
direction of matrix flow, and this is of particular significance at the 
seals where the direction of fluid flow reverses. 
Consider for the moment only the heat exchange element on the side 
of Cmax, Fig 2(a). For simple unmixed cross flow with uniform tempera- 
tures across the inlets and average temperatures at the outlets, the 
heat transfer rate for this element is 
Q = Cmax(Txi - Txo)(1/Nx) me a hs | 
Q = Cy(Tro - Tri)(1/N7) CIVEM V y 
Q = (hÀ)k ATayg(1/N HN) EN oues EE 
For & small enough element the arithmetic-mean temperature differ- 
ence may be assumed valid, so thet 
АЕ Tyo) (0/9 T. + Tro) > ee- [5] 
By using equations [2] , [5] 3 [u] , and E the outlet temperatures 
of the two streams may be solved for in terms of the inlet temperatures 
and put in the form (Appendix 1) 
Bu a Ne РАУ 
D MEET μι ο αμ. ο. 


By the same method the outlet temperatures of the element on the 


Tro 


side of Супр may be solved for in terms of its inlet temperatures giving 
Ino = Thi + Ка (Thy = Thi ) 9 ° € οϱ ο . ο 9 9 9 e e e (та | 
Tro = T; - Ka (Tri = Int ) 9 e е a 9 o © 9 o 9 e ° ° [το] 


By using the definitions of the parameters and proper algebraic 





manipulation the foregoing constants may be put directly in terms of 
the dimensionless parameters and the number of subdivisions of the matrix 


and fluid streams, giving 





= 
2м. (АА) 

К = 2 || Zee ms = 
| — кти, [+ WA] Ee й 
Ν n" mun 一 一 一 RE ле ὅ ο 

Ko 2, || Cus Cnn he M TUS [1 (hA)] 
2 М B 
E +c. м oo с ' 
з E з ura, [179 0] 
Cr | 
oF MA END C min 
= 2 ma ΠΤ, ΠΠ 
Κι 4 Cmin Nr NTU, [1+ (ва)! Шш 


Now consider the schematic representation of the composite heat 
exchanger made up of elements as shown in Fig. 3. For illustration the 
fluid streams and the matrix stream have each been divided into three 
equal substreams to form the elements. The double lines represent the 
area of the seals. It should be noted that the left edge is physically 
the same as the right edge and therefore the matrix inlet temperature of 
& particular matrix substream on the left must be identical to the matrix 
outlet temperature of that substream on the right. This will be referred 
to later as the reversal condition. 

In Fig. 3 where individual element temperatures are designated with 


the same label it is not implied that they are equal numerically, except 
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where indicated as constant across the fluid inlets. This is done to 
avoid additional subscript notation. Also, it is not to be implied from 
Fig. 3 that there is necessarily any symmetry with respect to the areas 
involved. 

Now for calculation purposes if the fluid inlet temperatures are 
given some convenient values, such as one and zero, and a matrix inlet 
temperature distribution is assumed on the left edge, the remaining 
temperatures of each element can be calculated in the order indicated 
by the element numbers by repetitive use of equations [6] and then [ τ/ è 
If the temperature distribution assumed on the left edge was correct 
then it would be duplicated on the right, i.e., the reversal condition 
was fulfilled, and the problem would be solved for the particular set 
of parameters used. If, however, this is not the case, then the result- 
ing temperature distribution on the right is now used on the left and the 
procedure repeated until the reversal condition is met. 

After several iterations of this type the cyclic character of the 
calculations for each element and column on a side becomes apparent and 
the beauty of automatic computing equipment is fully appreciated. It 
can be seen that if, in the process of computing the outlet temperatures 
of element one, the Ty, is put where the T4, was and a control number 
increased by one so that the T,; of element two.will be used next the 
same sequence of operations can now be used to compute element two, and 
so on down the first column of elements. Also, if the T,, of each ele- 
ment is put where ite То: was, then upon reaching the end of the column 
and “zeroing” the above control number the same cyclic operation can be 
used for the next colum. A second control number is used to indicate 


when one side is finished and then a similar procedure is followed on 
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the other side. In this scheme it is only necessary to provide stor- 
age for the Ty; and Tro of each row, which will be used for comparison 
purposes to determine if the reversal condition is met, and the Ty, of 
each column on one side and the Tr, of each column on the other side, 
which will be used to make a heat balance and compute the effective- 
ness when the reversal condition is met. 

That the reversal condition eventually will be met by this iterative 
scheme can be shown (Appendix 1) to depend on the physical conditions of 
the problem and the second law of thermodynamics. This will be called 
the condition for convergence and for this problem, using equal sub- 
divisions of each stream, it may be stated in the following form: 

NTU, [1+ 3/(nA)'] Cmin/Cmx £ 2 ....... [9] 

Accepting for the moment the idealizations and assumptions made, 
then the accuracy of the solution obtained, by the very nature of a 
finite-difference process, depends only on the number of subdivisions 
used. It is apparent from condition [9] that a greater number of sub- 
divisions will also enhance the convergence. 

A considerable number of subdivisions can be handled by automatic 
computing equipment to obtain a TORN accuracy and aid convergence, 
but the number of subdivisions used must be arrived at by a compromise 
between the accuracy desired and time available. The easiest way to 
illustrate the accuracy obtained for a given degree of subdivision is 
to solve the problem for several values of subdivision and extrapolate 
the results obtained to zero element area. As an example, consider the 


following results for the indicated parameter values: 


Cnin/Cmax = 0.90, NTU, = 6, ФЕ, in = (nA)! = 1,0 
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No. subdivisions Effectiveness 


8 0.7816 
16 0.7800 
32 0.7796 


These results are plotted versus element area and extrapolated to 
zero in Fig. 3a where it is found that the indicated effectiveness is 0.7795. 
This correction is not constant for all ranges of parameters, being 
largest for low values of C,/Cyin and high values of NTU, , the fore- 
going example being typical. The procedure used in this work was to 
calculate Tables 1 through 9 with 16 subdivisions, then recalculate 
Tables 1 and 9 with 32 subdivisions. Table 9 was calculated to show 
the effect of varying the parameter (hA)', and it also contains rep- 
resentative values from the other tables so that with these data all 
of the tabulated data were extrapolated to zero element area. 

The computer used has the equivalent of nine significant decimal 
figure accuracy; it was specified that the reversal condition be ful- 
filled to only five figures. Of more significance than this was the 
final requirement that before a solution was considerd as being obtained 
the heat balance error, which is a direct measure of error in the effec- 
tiveness introduced by accuracy of computation, be less than 0.024 per 
cent. 

The remainder of the computer program was essentially a system of 
"feeding" it the next set of parameters after each solution. Appendix 
2 describes the computer program and its various modifications in some 
detail. Several interesting features incorporated in it are: 

a. Up to 64 subdivisions could be accommodated without external 


memory storage simply by modifying certain control commands. 
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C. 


Final matrix temperature distributions were punched on cards 
for future use. 

A logical choice of initial estimate for matrix temperature 
distribution for each new problem was made from data available 
from work up to that point. 


Optional read out of all or selected temperatures was provided. 
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4. Comparison of Results 

The results obtained are reported in tabular (Tables 1 - 9) and 
graphical form (Figs. 4 - 11). Several general observations can be made 
from the data: 

a. As Cmin/Cmax decreases, for given values of the other parameters, 
the effectiveness increases. This is predicted by the limiting solution 
Lor O/C nin = 00 , equation TE and was pointed out 1n the illustrative 
problem of reference (2). This is illustrated in Fig. 12. 

b. By using NIU, as defined in the nomenclature, which contains 
(hA)' explicitly, the influence of (hA)' at a given NTU, is very small 
for Cmin/Cmax 2 0.9 . This conclusion was also reached by Coppage and 
London (2) in their evaluation of the work by Hausen (6), Saunders and 
Smoleniec (7), J.E. Johnson (8), and Iliffe (10). For Cmin/Cmax < 0.9, 
however, the influence of (hA)' becomes increasingly pronounced, and the 
variation for (hA)' = unity and 0.25 may be as mich as seven per cent for 
Cmin/Cmax = 0.1. See Table 9. 

ο. The difference in effectiveness for Cy/C,i, > 5 and the value 
given by the limiting solution for CL./Cgin = © is very small. This 
may be seen from inspection of Tables 1 through 8 or Figs. 4 through 11. 

R.W. Johnson (11) used much the same method employed here but had 
limited access to a computer and was therefore able to calculate only a 
relatively small number of the points reported here. Also due to limited 
computer time, at most ten апа in some cases only five subdivisions were 
used with no extrapolation to zero element area correction. Allowing for 
this the corresponding results are practically identical. 

Coppage and London (2) recommended the results of J.E. Johnson (9) 


which were obtained for Cmin/Cmax = unity. These results, in terms of 
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the parameters used here, were later incorporated in reference (12) with 
extrapolation to C un/Cmax = 0.7 . The variation in corresponding values 
of Table 1 of this work and Table 6, page II-32, of reference (12) which 
are for Cmin/Cmax = unity is less than one per cent. Only when Grn) Cras 
= 0.7 and Cr/Cmin = unity, where the extrapolation method used was ad- 
mitted to be most uncertain, is the variation slightly more than two per 
cent. 

An expression for effectiveness in closed form may be desir: able for 
certain uses such as computer analysis of plant performance. A relative- 
ly simple empirical formulation is suggested in reference (12) for 
Cmin/Cmax 2 0.7 and Cr/Cnin 2 2 and is of the form 


Eco - E AE 


E. E Eo Ks Πο τα د ر‎ O [10] 





where Kg = 1/9 ,m= -2, and En is determined from equation [1] . Such 
a form cannot be applicable for all values of parameters with single val- 
ues of Ks; and m, but it agrees reasonably well with the results reported 
here for the ranges specified. It has been found that with slight var- 
iation of K, and m very close agreement can be obtained in a small range 
of parameters which may be of particular interest. As an example, in 
the range 

πο ος. Z 0.90 


1.25 < pin < 5.0 


IN 


NIU. < 9.0 


5.0 o 


with Ks = 1/9 and m= -1.87, the error is less than one per cent. Fig. 


15 illustrates the method used to determine KS and m in this range. 
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4, Limitations 

The most questionable of the idealizations made is that of zero 
matrix conductivity, particularly in the direction of fluid flow. Α]- 
though the effect of matrix conductivity could be included in the finite- 
difference equations the additional degrees of freedom would make adequate 
coverage of all the parameters extremely voluminous and time consuming. 
If for a particular exchanger design the matrix conductivity was consider- 
ed important then the problem could be solved individually on a computer 
with adequate capacity. A simpler alternative is to neglect the conduct- 
ion and correct the effectiveness for this effect. From theoretical con- 
siderations Professor A.L. London has shown (13) that a correction factor 
for conductivity in the direction of fluid flow for approximately equal 


fluid capacity rates, of the form 


E M MM E ο ου 
E s 


should result in a somewhat pessimistic prediction of the reduction in 
effectiveness. (A in this case is the matrix cross-sectional area avail- 
able for heat conduction in the direction of fluid flow on both sides, 

and L is the flow length). Schultz (14) has solved for the effectiveness 
as a function of the conduction parameter in equation [u] for three cases. 


In the nomenclature used here, these three cases are for the folloving 


conditions: 
a. Cmin = Cmax » Cr/Cnin very large, k finite 
b: Cmin = Cmax ; Cr/CQ44 7 ?9 , k finite, regenerator sub- 
divided into sections in direction of fluid flow 
er Cmin = Cmax , Cr/Cmin finite, k= © , regenerator sub- 


Givided into sections in direction of fluid flow. 
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From the curves in reference (14) which allow direct comparison with 


the results presented here it is found that for 


КА Моо СЕТЕ 


a 


LC 
the ređuction in effectiveness is at most about one per cent. Addition- 
al work, not available at the time of writing, concerning the influence 
of conduction is reported in reference (15). | 
The effect of leakage can be calculated separately resulting in a 
correction to the effectiveness without leakage. This has been shown to 


be small for as mich ten per cent leakage (1). 
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5. Summary and Conclusions 
The numerical procedure used in this work assumes no symmetry, 
either with respect to heat capacity rates of the fluids or area ratios 
of the two sides, and therefore is not limited in this respect. The 
ranges of parameters which have been covered effectively are 
ο ο. > 0.10 


1.0 > (bA)' 0.25 
The results are tabulated at frequent enough intervals and over suffi- 
cient ranges of the parameters to include the gas-turbine regenerator 
appllcation as well as other applications, such as intercooling, where 
the values of the parameters might be quite different. 

The error introduced by the finite-difference type of solution has 
been minimized by use of sufficient subdivision and further extrapolation 
to zero element area so that four significant figures in the results are 
216164. 

The effects of deviations from the assumed idealizations have not 
been included in the results. 

Even though the results obtained cannot be made to fit a completely 
general closed form, a relatively simple expression for effectiveness 
can be obtained which will give good agreement over a reasonable range 


of parameters. 
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TABLE 1 


PERIODIC-FLOW HEAT 


Cain/Cmax 
(ВА)! 


EXCHANGER EFFECTIVENESS 


= 1.0 


Ξ 1.0 


END > —— Cr/ Chin Fa: 


1 

NIUo 

6.5 0.3221 
1.0 0.4665 
1.5 0.5478 
2.0 0.6007 
2.5 0.6385 
3,0 0.6672 
3.5 0,6900 
4.0 6 
4.5 0.7242 
5.0 0.7575 
5.5 0.7491 
6.0 0.7992 
6.5 0.7682 
m9 0.7165 
5 70.7855 
8.0 0.7901 
8.5 0.7961 
9.0 0.8017 
9.5 0.8068 
10.0 0.8115 


1 o 


0.4779 


0.6231 


0.6959 


0.7410 


0.7723 


0.7956 


0.8138 


0.8285 


0.8407 


0.8510 ˆ 


129 


0.3282 
0.4845 
0.5757 
0.6360 
0.6792 


0.7119 
0.7377 
0.7586 
0.7760 
0.7907 


0.8034 
0.8144 
0.8241 
0.8327 
0.8404 


0.8473 
0.8536 


0.8593 
0.8646 


0.8694 


2 


0.3304 
0.4912 
0.5861 
0.6491 
0.6941 


0.7280 
0.7546 
0.7760 
0.7957 
0.8086 


0.8213 
0.8323 
0.8419 
0.8504 
0.8579 


0.8647 
0.8708 


0.8763 
0.8814 


0.8860 


20 


0.4960 


0.6587 


0.7400 


0.7889 


0.8217 


0.8453 


0.8632 


0.8772 


0.8884 


0.8978 


BEES ri 
5 oo 

0.3329 0.3333 
0.4985 0.5000 
0.5978 0.6000 
0.6638 0.6667 
0.7110 0.7142 
0.7463 0.7500 
0.7739 ο με 
0.7959 0.8000 
0.8139 0.8181 
0.8290 0.8333 
0.8417 0.8461 
0.8526 0.8571 
0.8621 0.8667 
0.8704 0.8750 
0.8777 0.8823 
0.8842 0.8889 
0.8900 0.8947 
0.0995 0.9000 
0.9002 0.9047 
0.9045 0.9091 





10 


PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 


TABLE 2 
ج‎ 
1 1.25 

0.4714 0.4833 

0.6082 0.6319 

0.6760 0.7064 

0.7183 0.7528 

0.7477 0.7849 

0.7698 0.8090 

0.7872 0.8277 

0.8012 0.8429 

0.8129 0.8554 

0.8229 0.8661 





Cnin/ mex = 0.95 


(hA)' 


Cr/Cmin 


1.5 


0.4901 
0.6454 
0.7234 
0.7715 
0.8047 


0.8292 
0.8481 
0.8655 
0.8751 
0.8861 


- 1.0 


0.4970 
0.6592 
0.7406 
0.7903 
0.8240 


0.8484 
0.8672 
0.8820 
0.8940 


0.9040 


21 


0.5021 
0.6693 
0.7533 
0.8039 
0.8379 


0.8624 
0.8808 
0.8953 
0.9069 
0.9165 


0.5047 
0.6747 
0.7600 
0.8113 


0.8455 


0.8700 
0.8881 
0.9027 
0.9142 


0.9235 


0.5063 
0.6778 
0.7640 
0.8158 
0.8503 


0.8750 
0.8931 
0.9077 
0.9191 


0.9284 





NTU, 


10 


TABLE 3 


E Cr/Cmin 


1 


0.4763 
0.6156 
0.6846 
0.7275 
0.7513 


0.7795 


0.7970 
0.8112 


0.8230 
0,8331 


PERIODIC-FLOW HEAT 


1.25 


0.4887 
0.6405 
0.7167 
0.7642 
0.7970 


0.8216 


0.8407 
0.8561 


0.8689 
0.8797 


EXCHANGER EFFECTIVENESS 





Cnin/Cmax = 0.90 

(hnA)' » 1.0 
15 2 3 5 
0.4958 0.5030 0.5082 0.5110 
0.6547 0.6693 0.6800 0.6856 
0.7347 0.7530 0.7664 . 0.7755 
0.7841 0.800 0.8185 0.8265 
0.8181 0.8385 0.8534 0.8615 
0.8432 0.8636 0.8784 0.8865 
0.8625 0.8827 0.8972 0.9051 
0.8780 0.8978 0.9117 0.9195 
0.8906 0.9099 0.9234 0.9309 
0.9012 0.9199 0.9329 0.9401 


22 


0.5126 
0.6889 
0.7777 
0.8310 


0.8664 


0,8916 
0.9102 
0.9246 


0.9559 
0.9450 





TABLE 4 PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 


Cain Caas ~ 0.80 


(hA)' = 1.0 


六 


A. 1 1.05 1.5 2 2 oo 

0.1861 0.4994 0.5071 0.5149 0.5207 0.5251 
2 0.6299 0.6575 0.6750 0.6891 0.7011 0.7109 
3 0.7007 0.7364 0.7565 0.7769 0.7919 0.8043 
4 0.7443 0.7853 0.8077 0.8299 0.8461 0.8597 
5 0.7743 0.8190 0.8427 0.8655 0.8818 0.8957 
6 0.7967 0.8440 0.8683 0.8910 0.9069 0.9206 
7 0.8140 0.8634 0.8879 0.9101 0.9255 0.9386 
8 0.8279 0.8788 0.9035 0.9249 0.9392 0.9518 
9 0.8394 0,8914 0.9157 0.9565 0.9500 0.9619 
10 0.8492 0.9021 0.9260 0.9460 0.9586 0.9696 


25 





TABLE 5 PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 


Cmin/Cmax = 0.70 


(nÀ)' = 1,0 





ENS c ^ Cr/Cmin 


P. 1 1.25 5 2 2 за 

1 0.4959 0.5108 0.5187 0.5271 0.5555 0.5304 
2 0.6437 0.6737 0.6910 0.7087 0.7219 0.7326 
5 0.7155 0.7548 0.7770 0.7995 0.8160 0.8295 
4 0.7590 0.8044 0.8291 0.8535 0.8710 0.8855 
5 0.7887 0.8382 0.8643 0.8891 0.9063 0.9207 
6 0.8104 0.8629 0.8895 0.9139 0.9303 0.9439 
1 0.8272 0.8818 0.9085 0.9320 0.9473 0.9598 
8 0.8405 0.8967 0.9232 0.9456 0.9596 0.9709 
9 0.8514 0.9087 0.9348 0.9560 0.9687 0.9788 
10 0.8606 0.9188 0943 0.9642 0.9756 0.9845 


24 





NTU, 


M 


ο 0 - σι 


TABLE 6 


PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 


Cmin/Cmax Е 


(hÀ)' » 1.0 


0.50 


lE — е аи 


1 


0.5156 
0.6690 
0.7402 
0.7817 
0.8089 


0.8284 
0.8431 
0.8546 
0.8639 
0.8716 


1.25 


0.5322 
0.7045 
0.7869 
0.8555 
0.8670 


0.8894 
0.9062 
0.9191 
0.9293 
0.9376 


1.5 


0.5417 
0.7249 
0.8133 
0.8645 
0.8975 


0.9202 
0.9366 
0.9489 
0.9582 
0.9655 


2 


0.2215 
0.7461 
0.8399 
0.8928 
0.9256 


0.9469 
0.9615 
0.9717 
0.9789 
0.9842 


2) 


2 


0.5588 
0.7618 
0.8590 
0.9122 


0.9436 


0.9629 
0.9752 
0.9833 
0.9886 
0.9921 


QO 


0.5647 
0.7746 
0.8744 
0.9274 
0.9572 


0.9745 
0.9847 
0.9908 
0.9944 
0.9966 





NTUS 


Wi 


ο 0 - с 


TABLE 7  PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 
Cmin/Cmax = 0.30 
(hA)' = 1.0 
μα == Cr/ Cmin 
1 1.25 155 2 
0.5346 0.5558 0.5608 0.5761 
0.6902 0.7315 0.7556 0.780% 
0.7577 0.811% 0.8422 0.8729 
0.7951 0.8558 0.8895 0.9214 
0.8191 0.8838 0.9183 0.9490 
0.8560 0.9031 0.9374 0.9658 
0.8487 0.9171 0.9506 0.9764 
0.8588 0.9278 0.9602 0.9835 
0.8671 0.9363 0.9675 0.9881 
0.8740 0.9433 0.9729 0.9913 





26 


— — —— 
2 OO 
0.5846 0.5915 
0.7988 0.8136 
0.8945 0.9110 
0.9420 0.9566 
0.9670 0.9787 
0.9808 0.9895 
0.9885 0.9948 
0.9930 0.9974 
0.9955 0.9987 
0.9971 0.9994 





NTU, 


я + Ww N 


ο wo — ο 


TABLE 8 


PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 


Cri n/Cmax = 0.10 


(hA)' = 1.0 


一‏ و م 


1 


0.5534 
0.7070 
0.7680 
0.8009 


0.8223 


0.8378 
0.8498 
0.8595 
0.8675 
0.8743 


1.25 


0.5752 
0.7542 
0.8279 
0.8667 
0.8908 


0.9075 
0.9200 
0.9298 
0.9516 
0.9441 


1.5 


0.5877 
0.7820 
0.8630 
0.9042 
0.9283 


0.9440 
0.9550 
0.9632 
0.9694 
0.9743 


2 


0.6008 
0.8109 
0.8981 
0.9397 
0.9618 


0.9745 
0.9824 
0.9874 
0.9908 
0.9932 


2 


0.6104 
0.8319 
0.9219 
0.9615 
0.9802 


0.9892 
0.9939 
0.9964 
0.9977 
0.9984 


CO 


0.6186 
0.8487 
0.9391 
0.9753 
0.9900 


0.9959 
0.9983 
0.9993 
0.9997 
0.9999 





TABLE 9 


(БА)! 


1,00 


0.50 


. 00 


H- 


O 


«20 
125 


o 


. 00 


== 


H 


.00 


о 


.50 


1.00 


0.50 
0.25 


PERIODIC-FLOW HEAT EXCHANGER EFFECTIVENESS 


ο ο αρ ποτ Cr/Cmin ни! | 


1 


Cmin/Cmax 


0.6672 
0.6676 
0.6684 


0.7592 
0.7597 
0.7608 


0.8017 
0.8021 
0.8032 


Cnin/Cmax 


0.7007 
0.6960 
0.6919 


0.7967 
0.7902 
0.7848 


0.8394 
0.8323 
0.8263 


(ο 


2 


= 1.0 


0.7280 
0.7282 


0.7284 


0.8323 
0.8323 
0.8324 


0.8763 
0.8763 
0,8763 


= 0.80 


0.7769 
0.7752 
0.7742 


0.8910 
0,8888 
0.8872 


0.9365 
0.9344 
0.9328 
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1 2 
Cmin/Cmax = 0.95 
0.6760 0.7406 
0.6752 0.7404 
0.6751 0.7404 
0.7698 0.8481 
0.7687 0,8481 
0.7685 0.8478 
0.8129 0.8940 
0.8115 0.8936 
0.8111 0.8933 
Cmin/Cmax = 0.70 
0.7155 0.7995 
0.7075 0.7964 
0.7001 0.7943 
0.8104 0.9139 
0.8000 0.9103 
0.7908 0.9074 
0.8514 0.9560 
0.8402 0.9529 
0.8306 0.9503 


1 2 
Cmin/Cmax = 0.90 
0.6846 0.7530 
0.6826 0.752, 
0.6813 0.7521 
0.7795 0.8636 
0.7769 0.8628 
OTOL 0.8622 
0.8230 0.9099 
0,8198 0.9091 
0.8116 0.9084 
Cnin/Cmax = 0.50 
0.7402 0.8399 
0.7245 0.8331 
0.7099 0.8270 
0.8284 0.9469 
0.8102 0.9398 
0.7953 0.9334 
0.8639 0.9789 
0.8463 0.9735 
0.8327 0.9684 


TABLE 9 


(БА)! 


1.00 
0.50 


0.25 


1.00 
0.50 


0.25 


1.00 


0.50 
0.25 


CONTINUED 


Nw 7 cr/ Chin τμ - 


1 


2 


Cmin/Cmax = 0.30 


0.7577 
0.7338 
0.7131 


0.8360 
0.8130 


0.7957 


0.8671 
0.8472 
0.8327 


0.8729 
0.8608 
0.8484 


0.9658 
0.9546 
0.9439 


0.9881 
0.9808 


0.9754: 


1 


Cmin / Cmax 


0.7680 
0.7373 
0.7136 


0.8379 
0.8133 
0.7957 


0.8675 
0.8472 
0.8326 


2 


= 0.10 


0.8981 
0.8787 
0.8579 


0.9397 


0.9908 
0.9821 


0.9738 
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9. Appendix 1 - Derivation of Iterative Equations 


Equations [2] , D [4] and [5] eoncerning the element on the 


side of C ax, Fig. 2(a), from the body of the thesis are repeated here 
for convenience. 


το πο πο M .12) 
ο 
Q = Cy Ces - Tri) CLM) "um Lx 
Edu ο οι: Ὁ (4) 
AT, ο о) о (Та ть) ol [5] 


If ATavg from equation [ 5] is substituted in equation [ 4] and the 


value for Q thus obtained is substituted in equation [2] the following 
will be obtained: 


(A) (hab (Tri + The - Tez -Teo) (We Ix) = Cmax (Tat Т) СИ) 
This may be solved for Tyo, giving 


= 


БА), AN, Cmax ~ ChA) — Т 
ШЕ сет T E — iro И 
хо = [2M Cay + Coady J Te ом AN Coat Ах) 


In a similar manner T o may be solved for by using equation [5] , giving 


ç 


EN a 1, — [aM Ce c GA] s, s [It] 
Tro = | TREO, | o) E Су + (ha) й 


At this point four relations can be derived which will be useful 


in subsequent manipulations to get expressions in terms of the dimen- 


sionless parameters. From the definition of NTU, it is obvious that: 
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C min ; 


A 


An expression for reciprocal NTU, may be derived by multiplying and 


dividing it by NTU, as expressed in equation [15] , thus, 


ο μμ. ο ЕЕ САА), 
(MA), NTU, (hA)x NTU, [H8 (R4). — Cs 
ра)! : 


Ντι [1+ (nA) le, 


The other two.expressions, the derivations of which are obvious, are 


as follows: 








s ш κο [17] 
CM NTU, [1+ (nA)' ] 

and, 
Cr CHAN ew De 
(has — wr [rt CA) ] ο. 


To obtain ve in terms of only the inlet temperatures substitute 
equation [14] in equation [15] and solve for To. Starting with equa- 


tions [15] and [14] slightly rearranged, the algebra involved follows. 


T I = TT er 
fois tn] o t Ti T бт πι о [nd 
і б 2 = 
M | х Та) τὶ! - с [т oo, [da] 
X η 
(Αγ. 


Substituting [15a] in [13а] , 





| 2 
AN. ES | N 8 س‎ 6 і 
τα ΤΙ Ἴμα у | т ть 7 аст |+ а a 
ECT 十 MA (nA). 
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Each of the terms in the brackets may be put over the same common denom- 


inator which may then be multiplied out, giving 


«Cr 
4 Ne Na Ge „am „ 2Νκίς w z че. Nile ια ar. η нти 
утіх (нах | МУ (һа), Ko NTH (HA), Мт 办 


\ 4 Nu C+ | , 
十 “| 一 -一 一 一 2-2 17 
| (кА), Τ rc 


Multiplying through by (hA)'/2N,C, and adding and subtracting 2T,ı on 


the right gives 


AN: САХ Ne 20у 4 Cmax Ny Ny f | 
m CU M м |5, E [us C M. «ijr E Au "Ти ) 


Now solving for T,. gives a form convenient for iterative calculations, 





thus, 
— x Cmax Me AN, й 
D. = |х; -2)1* Ec NE 7 十 E Те m x M 


The term in brackets is constant for a given problem so that 


T yo B ο т К, ae Tos) | συν .: КОШИ: [ea] 
where K, may be put directly in terms of the dimensionless parameters 
giving 


г Y y "о 一 
| 4 Nc (5A J РЯ 


K = ¿| во ος Ps у ru] ον Eon 
er Cos Cmax 
To obtain T,, in térms of only the inlet temperatures substitute 
equation [13] in equation [14] and solve for T pọ- Starting with 
equations [13] and [1%] slightly rearranged, the algebra involved 


follows. 





N. й 2 2 
Е ey | (Ti + Ta) t n LN Τκι SLM | 36) 


NTUx 


(Ax, 


hh 





М, Се . =, ,- М Си і 
es 4] Tro 二 yi + lxo Т Е | e O, ΠΠ 
Substituting [15] in fw] , 


Qn, Cr | 7 SMC, | 
maana — + 一 = mi umge у - й 
| (АА); | ANY + | | E Шш | MA y + aM А | Hx 


NTUy 








Multiplying out the common denominator of all terms in brackets gives 


М Cr AN, 
4M Με Ον ДМС у any š _ üN. Mç G + «Να n 2 Jh | 
E. Ü (hA) * Wil. x (n AY УТИ» (на)х NTUx γι 





X 


+ [155 2-2) Txi 
NTU 


Multiplying through by NTU,/2N, and adding and subtracting 2T,, on the 


right gives 
[ 2W, Ce A πο , Μας + | T 
(hA) С гах ro (hA), Ny Cmax 


+ A (Tri πο 


Now solving for Tyo gives a form convenient for iterative calculations, 


thus, 


T Cy Νκ a Nx Cr 7-7] ~ 
o = ты жаи аад + о | Cue) Be] 


мах 


The term in brackets is constant for a given problem so that 


m p T Ka (TE πο | | | τ... К leb] 


νο, - 
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where Ko may be put directly in terms of the dimensionless parameters 





giving 
i cx E 
| ам, (па) 2 
М = 2 / i l min Cy Vy PE x cm | | | | [en] 
= max Coin ivt мт, [1+ (вА) J 


For the element on the side of C,;, , Fig 2(b), a set of equations 


similar to equations [e] througn [5] can be written. They are: 


ο ο Eee 
= NP 
доз ех (Ти = Tre) (Yale) 
ea] 
Q = (ha), May (Wr /ν.) "αν | 
-T o 8] 
Mar = (AN Tr: +T) E (tr) (Tri d Tao) | | і 
Following the same procedure as before results in two equations 
similar to equations [15] and [1] as follows 
= nA AN, m ES (hà), | [25 
Та ас ч πῆ. | (Tee +Tro) e } τ —— ^ 7 
and 
qM = ГОМ 0. = (кА) | [26] 
Tro Теа.) аж) 1 eet | PE 





These may be solved as before for the outlet temperatures in terms of 


only the inlet temperatures giving 


σα. --- ο. Nr ELS = i 1 
„та ој 2 m] Cu) oo τη 


lno 
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and 


& Ma y т. 
UR F Tv rafit smin Nr k | (Tai - ri) Ж . [23] 


The derivation of equations [27] and [28] is also readily appar- 
ent by comparison of equation [25] with 13) and [26] with [14] . The 
only difference is the substitution of the subscript "min" for "max" 
and "n" for "x" and then equation [27] follows from [19] and equation 

[28] follows from [20] . 


Equations [27] end [28] may be written as 





Tao = Tai + Ks (Tri = [m τ΄. ΠΣ [74] 
and 
Te = n EE Κι (Ту = E ͵ о 5 ΓΠΣ [76] 
where 
[ | ¿Nr =| 
ο EN вою ου τη 
Ка e us NTU, інф?) | ` 
and | 
aN, = -| 
Po Ne y min B. "ως: [$4] 
Κι Е | Cmin Ми мі Le 097] 


Equations [15] , [14] ; [25] and [26] are repeated here for con- 


venience. On the side of С пах : 


με. 


7 = _ ο Е 一 2M, Cmax - (ná msn 
[χο _ [= E (Tri * Γιο) y fcu + (ҺА) [т 1 b J 


(hA 
Tro т Кз ] г 2 № τ, za (hA), Ty! ei. 14 
E A | {πι}: Eu A 
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and on the side of Cmin : 


Τπο E “tn. | (Tvi í D y . EE (hA), Та Ki 


(ид), 


A A ΙΙ; E 
I" = ame + (oa), | Ti +) T PT ly vem, 


It is obvious that the first term on the right of each of the fore- 
going equations is a positive quantity. It can next be observed that 
the coefficient of the second term of each equation must also be positive 
for if it were not the outlet temperature in question would grow hotter 
as its inlet temperature became colder, or vice versa, which would be 
contrary to the second law of thermodynamics. 

This provides & convenient set of conditions to be met before the 
iterative scheme can be expected to converge. By using the relations of 


equations [15] through [18] these conditions may be stated as follows: 


NTU, [1+ Vinay pu < AN, JE) 
DEI a at i 
| ΝΙΝ Фо 


CHC min i 
Nu ji+ (nA | < 2⁄2 — s o5 [44] 


NTU] 1+ Cha) | 


ics Ж 
ο παν. 
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In the ranges of parameters considered and by using an equal num- 
ber of subdivisions of each stream it can be seen that condition [9] 18 
the most severe and therefore it alone can be used as the condition for 


convergence. 
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10. Appendix 2 - Computer Program Details 

The computer used for this work was the National Cash Register 102A 
general purpose digital computer and what follows assumes the reader has 
some familiarity with its operation. 

The basic arrangement of the computer program was designed to accom- 
modate any number of subdivisions of each of the three streams up to and 
including 64. It was found that, on the average, five iterations were 
required to meet the reversal condition and limit on the heat balance 
error with a judicious estimate of the initial matrix temperature dis- 
tribution. With 16 subdivisions of each stream the time required for 
each iteration was four minutes. For different degrees of subdivision 
the time required for each iteration would vary, as would be expected, 
as the square of the factor by which the subdivision was changed, e.g. 
32 subdivisions of each stream would require 16 minutes for each iter- 
ation. It was determined that 16 subdivisions of each stream could be 
expected to give accuracy to ο. five units in the fourth place of the 
effectiveness. To finish in a reasonable time and obtain the desired 
four place accuracy it was decided that all the results would be cal- 
culated using 16 subdivisions of each stream to be followed by suffic- 
lent check calculations with 32 subdivisions to permit extrapolation of 
the results to infinite subdivision. 


The use of the various cells of the memory was as follows: 


Cells Use 
0000-0477 Program 
0500-0577 Initial T's 
0600-0777 - Program 


1000-1077 Current parameters and K's, and 


temporary storage. 
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1100-1177 Outlet temperatures on side of Cmax 
1300-1377 Outlet temperatures on side of Chin 
1500-1577 Initial-Intermediate-Final T's 

The composite computer program flow diagram is shown in Fig. 14 and 
the contents of cells 0000-0477 and 0600-0777 are recorded in Table 10. 

Box 10 and 11 of Fig. 14 accomplish the iteration process. This is 
done by the repetitive application of equations [6] to each element on 
the side of Cmax in box 10 and equations [1] to each element on the side 
Of Cmin in box 11. Since the internal temperatures are of no concern 
in this problem it was only necessary to provide permanent storage for 
the fluid outlet temperatures of each column and the initial and final 
rotor temperatures of each row. 

Box 12 determines if the reversal condition is met by taking the 
difference between the initial and final T, of each row. A difference 
greater than the specified limit will cause the iteration cycle to 
start over. Also the largest difference is selected so that if the re- 
versal condition is met a means of manual inspection is provided, if 
desired, by causing this difference to be printed out upon putting switch 
2020 in the up position (box 17). 

If the reversal condition is met then the heat balance error is 
computed in box 16 and if this is greater than the specified limit the 
iteration cycle is started over. The equation used to compute the 


heat balance error was arrived at as follows: 


Qx 7 Qa 
Qx 





оо ы a = 


ο ae u Ti ) 
a (7; 一 ο 
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where Tno and Txo are the average fluid outlet temperatures. Since 
Ty, was taken as unity and T,; as zero then the heat balance error re- 


duces to 








Crain (т, ) 
error = / = Ema ο.” 
` / з Ee 
The average fluid outlet temperature on each side was computed by sum- 
ming the individual outlet temperatures of each column and dividing by 
the number of subdivisions. 
Another means of optional manual inspection was provided in box 19. 
If desired, the heat balance error would be converted from octal to decimal 
and printed out by putting switch 2040 in the up position. 
By definition, the effectiveness may be computed by either of the 
following equations: 


EE = у ) ει ] 
= ΕΕ TE TEC a 


Cmin Thi = Tex ) 


"DS Cc τη, = P SEEE E: [sib] 
m (πι == πι) 


where Tho aná To, are the average fluid outlet temperatures on the hot 
and cold sides, respectively. From the conditions assumed for calcula- 
tion purposes the following may be observed: 

Toi = Thi = O 

Thi = Txt = 1 

Co = Cmin , and, 

Ico = Tno 
and therefore the expression for effectiveness from equation [515] ге- 


duces to 


Е = (Tno)avg S. ...... .. 52) 


De 





(Tho ave will have already been computed in the process of calculating 
the heat balance error and therefore in box 21 it is only necessary to 
convert this to decimal and print it out. 

The final T,'s punched on cards in box 22 were used as estimates 
for the initial Ty's in subsequent calculations. This was usually done 
by using the final Ty's obtained at a given Cmin/Cmax as the initial 
T,'s for the next lower value of Cyin/Cmax, at corresponding values of 
the other parameters. 

The remainder of the program consists of a provision for interupting 
the calculations at logical stopping places and a system of inserting 
new values of parameters after each solution. These are self-evident 
in Fig. 14. Parameters are entered as floating point numbers and the 
arithmetic involved in calculating the K's is done in floating point. 
The K's, scaled down by a factor of two, are then converted to single 
precision and the subsequent iteration is accomplished in fixed point, 
single precision. 

The basic computer program which is recorded in Table 10 is set 
up to use 16 subdivisions of each stream. To modify it for 8 or 32 


subdivisions the following cells should be changed as indicated: 


Cell 8 Subdivisions 32 Subdivisions 
0307 21050021000010 2105002100004 0 
0403 5 了 
0406 35051021001510 3505402100150 
O410 36200015102001 36200015402001 
0411 35200021001110 55200021001 10 
О%2% 36150720002001 36153720002001 
0425 36150720031507 36153720031537 
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0430 35200021001310 35200021001 340 
0433 2ff0005 2ff0005 

0436 301310043 32007 30134004 332007 
0437 301110043 32007 401140043 32007 
0443 35151021000510 35154021000540 
0447 36051015102000 36054015402000 
0456 Ol, 3500030000510 Ol 300030000540 
0462 05300030001510 05300030001540 


In its present form the program by-passes box 25 and goes directly 


to box 27. 


get results with (hA)' = 1, 1/2, and 1/4 enter the following changes: 


0007 21100604020002 
0367 33100606760371 
0370 34300021000373 
0371 36100604001006 
0572 314300021000006 


In its present form the program produces results with C,/Cmyp = 


1, 1.25, 1.5, 2, 3, and 5. 


changes: 
0375 34067210020600 
0376 34067310030600 


This causes all calculations to be made with (hA)' 


To eliminate C,/Cmin = 5 enter the following 


Space is provided in the iteration subroutines for entering appro- 
priate commands to get any desired temperatures printed out before they 
are lost. These may be printed in decimal simply by placing in cell 
2000 and unconditional transfer to the decimal conversion and print sub- 


routine at cell 0700, 




















= — 
"чш ца 
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. 
А. 


FIRST 
COMMAND 
0000 


Set current Cnin/Cmax = 1 


0010 Compute K's x 
Test switch 2100 


ECAM Fill initial Tr's from cards 
0150 Initial Tr's double stored 
Iteration on slide of Cmax ` Down 


1 ЕЗИ Iteration on side of Cmin 
12 0314 Reversal condition met? No 


10 
1 


Yes 
Heat balance error small enough? 
Yes 


MO Οι 


S 
а 
Š 5 
+ 
Ὁ 


Test switch 2040 


0347 Convert E to decimal and print 
Punch final Tr'!s on cards 


Test switch 2010 Up 


š 


й 
ча 
мл 
го 


Down 


Set current Cr/Cmin = 1 Put in next Cmin/Cmax 


ram Set current NTUo - 1 30 | обоо Put in next Cr/Omin 
Set current (ha)! = 1 56 [овоз [Put in next кто 
ru Print current parameters Е 0571 |Put in next (ҺА)! 


Compute 


NO 


No 


NO 


No 


Í uu 


= | | jon rear 
(nA)! finished? 


` 


15 0306 | Print final Tr's and stop 


Up 


No 


17 0326 Test switch 2020 U Print lergest reversal 
0526 | а 18 | 0221 temperature difference 
Print heat balance error 
BEE 


Fig. 14 COMPUTER PROGRAM FLOW DIAGRAM 
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TABLE 10 


0000 31040004011000 
0001 31040004011002 
0002 31040004011004 
0003 31040004011006 
0004 21100004020002 
0005 21100204020002 
0006 21100404020002 
0007 34300021000010 
0010 31040004012000 
0011 31100610072002 
0012 34300021000733 
0013 31200420052000 
0014 31100410052002 
0015 34300021000725 
0016 31200420052002 
0017 31040304012000 
0020 34300021000760 
0021 31200420051040 
0022 31200420052000 
0023 31100610072002 
0024 34300021000725 
0025 31200420051042 
0026 31100210032002 
0027 34300021000725 
0050 31200420051044 


0031 31040004012000 


COMPUTER PROGRAM COMMANDS 


56. 


0032 34300021000760 
0033 31200420051046 
0034 31104210432000 
0035 34300021000725 
0036 31200420051050 
0037 31100010012002 
0040 34300021000760 
0041 31200420051052 
0042 31100210032000 
0043 34300021000725 
O0hh 3120042005105} 
0045 31200420052002 
0046 31040004012000 
0O47 34300021000760 
0050 31200420051056 
0051 31105210532000 
0052 31105610572002 
0053 343000210007 33 
0054 31200420051060 
0055 31200420052002 
0056 31040004012000 
0057 34300021000733 
0060 31200420052002 
0061 34300021000760 
0062 31200420051010 


0063 31105010512000 





TABLE 10 CONTINUED 


0064 31105410552002 
0065 343000210007 33 
0066 31200420052002 
0067 31040004012000 
0070 34300021000733 
0071 31200420052002 
0072 34300021000760 
0073 31200420051012 
007% 31104010412000 
0075 31104610472002 
0076 34300021000733 
0077 31200420052002 
0100 31040004012000 
0101 34300021000733 
0102 31200420052002 
0103 34300021000760 
0104 31200420051014 
0105 31104410452000 
0106 31100210032002 
0107 34300021000733 
0110 31200420052002 
0111 31040004012000 
0112 34300021000733 
0113 31200420052002 
0114 34300021000760 


0115 312004 20051016 


2T 


0116 
0117 
0120 
0121 
0122 
0125 
0124 
0125 
0126 
0127 
0130 
0131 
0132 
0133 
0134 
0135 
0136 
0137 
0140 
O141 
0142 


0143 


30101110101011 
30101310121013 
30101510141015 
30101710161017 
33210010110010 
33210010130010 
33210010150010 
33210010170010 
17210030000640 
343000210001 30 

35040421000131 
35052021001520 
35013104050131 
34040601310131 
35041221000157 
35041321000141 
35041421000146 


34300021000140 


35040121002000 
36200015202001 
30200104002001 


25200110112002 


0144 25200110132003 
0145 36200020022000 
0146 35152020031520 


0147 3430002100015} 





TABLE 10 CONTINUED 


0150 00000000000000 
0151 00000000000000 
0152 00000000000000 
0153 00000000000000 
0154 3501104070141 
0155 35014604050146 
0156 34041001410141 
0157 35200021001120 
0160 34300021000170 
0161 00000000000000 
0162 00000000000000 
0163 00000000000000 
0164 00000000000000 
0165 00000000000000 
0166 00000000000000 
0167 00000000000000 
0170 35015704000157 
0171 34041101570135 
0172 34300021000200 
0173 00000000000000 
0174 00000000000000 
0175 00000000000000 
0176 00000000000000 
0177 00000000000000 
0200 35042721000221 


0201 3504242100020} 


0202 35042521000211 
0203 35210021002000 
0204 36147720002001 
0205 30200104002001 
0206 25200110152002 
0207 25200110172003 
0210 35200020022000 
0211 36147720031477 
0212 34300021000216 
0213 00000000000000 
0214 00000000000000 
0215 00000000000000 
0216 36020404150204 
0217 36021104050211 
0220 34020404260204 
0221 35200021001320 
0222 34300021000231 
0223 00000000000000 
0224 00000000000000 
0225 00000000000000 
0226 00000000000000 
0227 00000000000000 
0230 00000000000000 
0231 35022104000221 
0252 5035002210201 


0233 34300021000314 





TABLE 10 CONTINUED 


0234 00000000000000 
0235 00000000000000 
0236 00000000000000 
0237 00000000000000 
0240 3504342100024 3 
0241 35043521000255 
0242 35210021002000 
0243 30132004332007 
O244 35200020072000 
0245 35024304150243 
0246 34043602430243 
0247 30200004001040 
0250 31210020002000 
0251 31100010012002 
0252 34300021000725 
0253 31200420052000 
0254 35210021002003 
0255 30112004332007 
0256 35200320072003 
0257 35025504150255 
0260 34043702550255 
0261 36040120032003 
0262 31210020032002 
0263 35200204002002 
0264 34300021000265 


0265 34300021000760 


20 


0266 30200520042005 
0267 36040120052005 
0270 30200504002005 
0271 34200504400300 
0272 34300021000326 
0273 00000000000000 
0274 00000000000000 
0275 00000000000000 
0276 00000000000000 
0277 00000000000000 
0300 35044221000300 
0301 35150021000500 
0302 35030104050301 
0303 34044303010301 
0304 17201030000306 
0305 343000210001 34 
0306 21044404450001 
0307 21050021000020 
0310 22000000000000 
0311 34300021000134 
0312 00000000000000 
0313 00000000000000 
O314 35210021001050 
0315 35044621000316 
0316 36052015202000 


0317 34200010500324 





TABLE 10 CONTINUED 


0320 35031604500316 
0321 34044703160316 
0322 34105004410300 
0%23 34300021000240 
0324 35200021001050 
0325 34300021000320 
0326 17202030000331 
0327 17204030000333 
0330 34300021000340 
0331 21105021000001 
0332 34300021000327 
0333 35200521002000 
0334 34300021000700 
0335 34300021000340 
0336 00000000000000 
0337 00000000000000 
0340 34300021000347 
0341 00000000000000 
0342 00000000000000 
0343 00000000000000 
0344 00000000000000 
0345 00000000000000 
0346 00000000000000 
0347 35104021002000 
0350 34300021000700 


0351 210444044 50001 


0352 35045721000353 
0353 05300030001520 
0354 13046004550000 
0355 30300004611077 
0356 13046004550004 
0357 30300004611077 
0360 35035304530353 
0361 34046203530353 
0362 34300021000363 
0363 34300021000364 
0364 17201030000366 
0365 34300021000367 
0366 22000000000000 
0367 34300021000373 
0370 00000000000000 
0371 00000000000000 
0372 00000000000000 
0373 34065210040603 
0374 34065310050603 
0375 34067410020600 
0376 34067510030600 
0377 34300021000611 
0400 00000000000000 
0401 004.00000000000 
0402 02000000000000 


0403 00000000000006 





TABLE 10 CONTINUED 


O404 35050021001500 
0405 00000100000001 
0406 35052021001520 
0407 00000000010000 
010 36200015202001 
0411 35200021001120 
0412 35200021001100 
0413 36200015002001 
0414 35150020051500 
0415 00000100000000 
0416 30050004002000 
0417 30052004002000 
0420 30152004002000 
0421 5015000002000 
0422 30110004002000 
0423 30112004002000 
O24 36151720002001 
0425 36151720031517 
0426 36147720002001 
0427 35200021001300 
0430 35200021001320 
0h31 30130004002000 
0432 30132004002000 
0433 02000000000004 
0434 30130004332007 
0435 30110004332007 
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0436 30132004332007 
O437 301120043 32007 
O440 00000100000000 
04%] 00000004000000 
0442 35150021000500 
O443 35152021000520 
OUL} 00341212341212 
01+5 10000000000000 
04146 36050015002000 
0447 36052015202000 
0450 00000100010000 
0451 04300030000500 
0452 15555555555556 
0453 00000000000010 
0454 00000000000500 
0455 11111111111111 
0456 04300030000520 
0457 05300030001500 
0460 26666666666666 
O61 00000000001000 
0462 05300030001520 
0463 00000200020000 
0464 31066406651002 
0465 31210006551000 
0466 00000000000000 


0467 00000000000000 





TABLE 10 CONTINUED 


0470 00000000000000 
0471 00000000000000 
0472 00000000000000 
0473 00000000000000 
0474 00000000000000 
0475 00000000000000 
0476 00000000000000 
0477 00000000000000 
0600 31066406651002 
0601 35060004630600 
0602 34300021000002 
0603 31040004012000 
0604 31100410052002 
0605 34300021000733 
0606 3120042005100} 
0607 31040004011006 
0610 34300021000006 
0611 35016421000600 
0612 31100010012000 
0613 35210021002002 
0614 36210006632003 
0615 34300021000733 
0616 34200521000621 
0617 35046521000621 
0620 22000000000000 


0621 31210006551000 
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0622 35062104070621 
0623 34300021000001 
0624 00000000000000 
0625 00000000000000 
0626 00000000000000 
0627 00000000000000 
0630 00000000000000 
0631 00000000000000 
0632 00000000000000 
0633 00000000000000 
0634 00000000000000 
0635 00000000000000 
0636 00000000000000 
0637 00000000000000 
0640 35045121000646 
0641 05300030002100 
0642 06045204550000 
0643 30300004541077 
0644 06045204550004 
0645 30300004541077 
0646 0300030000520 
0647 35064604530646 
0650 34045606460641 
0651 34300021000130 
0652 00000000000004 


0653 00500000000000 





TABLE 10 CONTINUED 


0654 00000000000001 
0655 00746314631463 
0656 00714631463146 
0657 00631463146315 
0660 00546314631463 
0661 00400000000000 
0662 00231463146315 
0663 00063146314632 
0664 00000000000001 
0665 00500000000000 
0666 00000000000001 
0667 00600000000000 
0670 00000000000002 
0671 00400000000000 
0672 00000000000002 
0673 00600000000000 
0674 00000000000003 
0675 00500000000000 
0676 02000000000001 
0677 00400000000000 
0700 27300007162001 
0701 35071420010713 
0702 35200021000722 
0703 05300030000714 
0704 32200207002000 


0705 2520022007 2002 
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0706 32200220012000 
0707 32210020012002 
0710 27200120052001 
0711 34200121000705 
0712 21200007200001 
0713 34300021000351 
0714 34300021000000 
0715 02000000000004 
0716 02000000000030 
0717 00500000000000 
0720 01000000000000 
0721 0070000000000 
0722 00463511721510 
0723 02000000000001 
0724 00000000000001 
0725 30300007162004 
0726 350714200407 352 
0727 35200020022004 
0730 25200120032005 
0731 3120042005200) 
0732 34300021000253 
0733 30300007162004 
0734 35071420040745 
0735 33200020020746 
0736 36200020022004 
0737 30200120042005 





- TABLE 10 CONTINUED 


O740 35200321002006 
O741 3520022100200 
0742 35200520062005 
0743 37200521000753 
0744 31200420052004 
0745 34300021000113 
0746 36200220002004 
0747 30200320042005 
0750 35200121002006 
0751 35200021002004 
0752 34300021000742 
0753 27200507232005 
0754 30200507232005 
0755 27200507242005 
0756 35200407242004 
0757 34300021000744 
0760 30300007162004 
0761 35071420040766 
0762 3620002002200} 
0763 23200120032005 
0764 37200521000767 
0765 31200420052004 
0766 34300021000266 
0767 27200507232005 
0770 30200507232005 
0771 27200507242005 


0772 35200407242004 
0773 34300021000765 
0774 00000000000000 
0775 00000000000000 
0776 00000000000000 


0777 00000000000000 
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